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Abstract 

The production rates of the orbitally excited D** mesons, and D*^, are measured 
with the 4.1 million hadronic Z decays recorded by the ALEPH detector during 
1991-1995. The D** mesons are reconstructed in the decay modes 

^ D*°K+ and ^ D°K+. The production rate of the is measured to be 
/(Z —> D^) = (0.52 ± 0.09 ± 0.06)%, under the assumption that the two considered 
decay modes of the saturate the branching ratio. The production rate of the 
is determined to be /(Z —> D*^) = (0.83 ± 0.29lQ'j’3)%, assuming that the 
branching fraction of the decay is 45%. The production rates in 

Z —> cc and Z —> bb decays are measured separately. 
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1 Introduction 


Four orbitally-excited D** mesons, with angular momentum F = 1, are expected to exist in 
addition to the pseudoscalar and vector mesons Dg and D*. In the Heavy Quark Effective 
Theory ii the spin of the light quark couples with the orbital angular momentum to 
give j = I or j = |. When coupled to the spin of the heavy quark, two doublets are 
obtained with the quantum numbers = O’*", and respectively. These 

states are expected to decay mainly into D K or D* K modes. 

Only certain S- and D-wave decays are allowed by spin and parity conservation. The 
states of the (j = |) doublet decay via S-wave transitions and are therefore expected to 
be broad (T > 100 MeV); they have not been observed yet. The (j = |) doublet states are 
much narrower since they can only decay via D-wave modes. The two narrow states^ 
and Dg^ have been observed by ARGUS and CLEO P, At LEP, mesons have 
been observed by OPAL [^]. The properties of the four D** states are listed in Table [^. 
The actual physical particles could be superpositions of the individual states. 

In this analysis the production rates of the two narrow states are measured in the decay 
modesQ —> D*+K°, ^ D*°K+ and —> D°K+. In Z —cc events, D** mesons 

are produced in the fragmentation of primary c-quarks, whereas in Z —>■ bb events they 
can only be produced in decays of b hadrons. To study these two contributions separately, 
b- and c-quark-enriched event samples are also used. 
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Table 1: Masses and decay modes of D** mesons. 


2 The ALEPH detector 

The ALEPH detector and its performance are described in detail elsewhere P |^. Only a 
brief overview of the apparatus is given here. Surrounding the beam pipe, a high resolution 
vertex detector (VDET) consists of two layers of double-sided silicon microstrip detectors, 
positioned at average radii of 6.5 cm and 11.3 cm, and covering respectively 85% and 
69% of the solid angle. The spatial resolution for the rep and z projections (transverse 
to and along the beam axis, respectively) is 12 pm at normal incidence. The vertex 
detector is surrounded by a drift chamber with eight coaxial wire layers with an outer 
radius of 26 cm and by a time projection chamber (TPC) that measures up to 21 three- 
dimensional points per track at radii between 30 cm and 180 cm. These detectors are 
immersed in an axial magnetic field of 1.5 T provided by a superconducting solenoidal 
coil and together measure the transverse momenta of charged particles with a resolution 
(j(pt)/pt = 6 X 10“^ pt © 0.005 (px in GeV/c). The resolution of the three-dimensional 
impact parameter for tracks having information from all tracking detectors and two VDET 

^The notation used follows that of the Particle Data Group |^] for the D** mesons. 

^Throughout this letter the notation used for particles implies the charge conjugate modes as well. 






























hits (a VDET “hit” being dehned as a space point reconstructed from the rep and z 
coordinates) can be parametrized as a = 25pm + 95pm/p (p in GeV/c). The TPC also 
provides up to 338 measurements of the specihe ionization of a charged particle. The TPC 
is surrounded by a lead/proportional-chamber electromagnetic calorimeter segmented 
into 0.9° X 0.9° projective towers and read out in three sections in depth, with energy 
resolution a{E)/E = 0.18/^/E + 0.009 {E in GeV). The iron return yoke of the magnet 
is instrumented with streamer tubes to form a hadron calorimeter, with a thickness of over 
7 interaction lengths, and is surrounded by two layers of muon chambers. An algorithm 
combines all these measurements to provide a determination of the energy flow |]^ with 
an uncertainty on the measurable total energy of <j{E) = + 0.6) GeV. 


2.1 Event and track selection 

This analysis uses 4.1 million hadronic events recorded by the ALEPH detector at centre- 
of-mass energies close to the Z mass. The events are selected with the charged particle 
requirements described in [^. The helix £t of the charged particle tracks used for the 
D** reconstruction must have a per degree of freedom smaller than 5, and their polar 
angle must satisfy | cos^l < 0.95. Finally the distance to the primary vertex in the plane 
transverse to the beam axis has to be less than 2 cm and in the beam direction less than 
10 cm. 


2.2 Particle identification 

Gharged kaons are identihed by means of the specihe ionization loss dE/dx in the TPG. 
The TPG provides two different measurements of the deposited energy, from the wire 
and the pad readout. The wire readout is only used if at least 50 individual wire 
samples are available. In this case the pad readout is ignored. If the wire information 
is insufficient the pad readout is exploited. In both cases the particle identification is 
based on the dE/dx estimator (?"k), defined as the difference between the expected 
and the measured ionization loss expressed in terms of standard deviations for the pion 
(kaon) mass hypothesis. A track is accepted as a kaon if its momentum is greater than 
1.5 GeV/c, Irxl is less than 2.5 and it satisfies < 0. 

Neutral kaons are reconstructed in the deay mode Kg 


TT ' TT 


as described in IIT2| . For 


the identification of the Kg their long lifetime (cr = 2.7cm) is exploited by only accepting 
neutral vertices with a distance of at least 1.5 cm to the primary vertex and a < 13 of 
the vertex fit. The reconstructed mass of the Kg candidate has to be within 12 MeV/c^ 
of the nominal Kg mass. 

Gharged pions are selected by requiring \r/ to be smaller than 2.5. Neutral pions are 
reconstructed by a mass-constrained fit to the detected photons in the electromagnetic 
calorimeter []T|. The momentum of the candidate has to be greater than 1.5 GeV/c 
and the of the mass fit less than 8. 


2.3 Selection of c- and b-quark-enriched event samples 

To determine the production rates of D** in Z —cc and Z ^ bb events separately, c- 
and b-quark-enriched event samples have been selected. Due to the good spatial resolution 
of the vertex detector, the probability ("Pevt) that all charged tracks originate from the 
primary vertex, determined with the algorithm described in , can be used to separate 
long-lived b hadrons (low probability Pevt) from c and uds events. 




















An additional separation between c and b events is obtained by a cut on the scaled 
energy xe (normalized to the beam energy) of the D** candidate, because less energy is 
available for c quarks originating from a b-hadron decay than in direct Z —»• cc decays. 

A D° is reconstructed in all the decay modes considered. In order to further separate 
light (uds) from heavy (c and b) quark events, a cut on the proper decay time of the 
candidate is performed. The proper decay time is determined from the decay length 
/, calculated as the distance of the D° decay vertex from the primary vertex, projected 
along the momentum vector of the D°, and taking t = M^o ■ l/p. 

The cuts used to select the enriched samples are listed in Table About 4% of all 
events are rejected because they pass the cuts for the c- as well as the b-quark selection. 


c-quark-enriched sample 

b-quark-enriched sample 

0.1 ps < t <1.5 ps 

t > 0.4 ps 

- logio(^evt) < 4 

-logio(^evt) >1-5 

Xe > 0.4 

0.25 <xe< 0.65 


Table 2: Cuts for the selection of the c- and b-quark-enriched event sample. 


3 Reconstruction of and D*^ 

3.1 Decay channel ^ 

The D*"*" mesons are reconstructed in the decay channel and D° mesons 

in the four decay channels D° —K“7r+, D° —K tt+tt", D° ^ K“7r’''7r“7r’'' and 
D° —> Because the combinatorial background is suppressed by the good 

resolution on the mass-difference between the and D° meson, the cuts for the selection 
of D° mesons are rather loose. 

To avoid double counting in the D° decay channels with charged kaons and pions, 
the track with the smallest value of tk is assumed to be the kaon track. A £t to a 
common decay vertex is performed for the tracks of the selected charged particles and 
the has to be less than 20. The momentum of the D° candidate has to be greater than 
7 GeV/c. Additional cuts are listed in Table ^ The selected D° candidates are combined 
with all charged pions with a momentum less than 3 GeV/c. For the D*^ candidates with 
a momentum greater than 8 GeV/c, the mass difference AMi = —m(D°) is 

calculated. The AMi distributions for the four D*’ decay channels are shown in Fig. [^. 
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Table 3: Guts for the reconstruction of D° mesons in the decay channels D° —> K vr’*', 
D° —>■ > K“7r+7r“7r''' and 




















Figure 1: Reconstructed mass-difference distribution of AMi. The quoted numbers of 
D*"*" candidates are obtained from a fit to the distribution. 

The number of reconstructed D*"*" mesons is determined by a fit to the measured AMi 
distributions. The signal is parametrized as a Breit-Wigner and the background shape 
with the following function: 


dN 

d{AM) 


oc (AM 


af ■e-c^^M-a)_ 


Each D*+ candidate inside the signal region (±1.5 MeV/c^) is combined with all Kg 
candidates with a momentum greater than 2.5 GeV/c. The combination of the D*"*" and 
the Kg is required to have a scaled energy xe greater than 0.25. The mass difference 
AM 2 = m(D*+,K°)—m(D*+)—m(K°)±mtrue(K°) is calculated, where mtrue(K°) is the 
nominal K*^ mass, because this quantity has the best resolution. The resulting mass 
difference distribution AM 2 is shown in Fig. The distribution shows an excess of 
events in the signal region around 525 MeV/c^ . The combinatorial background (Fig. |]b) 
is determined with the help of candidates from sidebands above the D*"*" mass window 
(Fig. i- The number of signal events and the mass difference AM 2 are determined by 
a simultaneous unbinned log-likelihood fit to the signal and background distributions, 
with the decay width of the fixed to 1.5 MeV. The signal is fitted with a Breit- 
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Figure 2: Reconstructed mass-difference distribution for the decays —>■ D*'''K‘’(a and 

b), D*°K+(c and d) and D°K+(e and f). The full curve is the result of an 

unbinned log-likelihood £t and the dashed curve represents the combinatorial background. 


Wigner function, convoluted with the detector resolution. The detector resolutions for 
the reconstruction of mesons in the different decay channels are determined 
from simulated events. The resolution is approximately 1.7 MeV/c^for the D° decay 
channels ^ ^ K^vr+Tr^and D° —> K“7r’''7r“7r’''and 2.6 MeV/c^ for the channel 

D° ^ K-tt+ttF 

The results of the £t to the AM 2 distribution are listed in Table The systematic 
errors on the mass differences are obtained by varying, in the £t, the assumed natural 
width of the meson and the detector resolution, as described in Section 

3.2 Decay channels ^ and D*^ —^ 

The D*° mesons decay into a D° by emitting a or a photon. Because of the higher 
background in these channels the selection of candidates has to be more restrictive 
than for candidates used in the reconstruction of mesons. An acceptable signal- 


5.1 
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m(K' 71 ^) [GeV/c^] 

Figure 3: Invariant mass of the D° ^ 'K~ candidates. The events in the shaded region 
are used for the reconstruction of the and D*^ ^ decays. 


Decay channel 

No. of events 

AM [MeV/c^] 

D+i ^ D*+K0 

40,9 ± 

525.3 ± O.Ostat ± O.lsyst 

D+3^D*°K+ 

51.3 ± 

528.7 ±1.9,tat±0.5,y,t 

d: 2 +^dok+ 

63.6 ± 21.9.,., 

704 ± 3stat ± Isyst 


Table 4: Results of the fits to the flavour-independent sample. 


to-background ratio is only achieved in the decay channel D° — Yr and therefore only 
this channel was used for the analysis. The D° candidates are rejected if the invariant mass 
of the D° and any 7 r+ in the event is within 2 MeV/c^of the D*"*" mass window. The cuts 
on the momentum of the candidate and on that of the decay particles are tightened to 
8 GeV/c and 3 GeV/c, respectively. The D° candidate is also required to have a proper 
decay time greater than 0.2 ps. The latter cut removes combinatorial background, peaked 
at small proper time. The invariant mass distribution of the selected D° candidates is 
shown in Fig. |^. The or photon from the decay is not reconstructed in this analysis. 

Instead the combination of the D° with all charged kaons in the event with momentum 
greater than 3 GeV/c is taken and the mass difference AM 3 = m(D°,K’'')—m(D‘^) is 
calculated. Because of the small Q-value in the decay of only 7 MeV, 

this procedure does not significantly change the resolution, which is determined from 
the simulation to be (2.9 ± 0.1) MeV/c^. In the case where the emits a photon the 
resolution is (7.1 ± 0.2) MeV/c^ . 

The meson is allowed to decay directly into the D° K"*" final state and 

therefore the mass difference AM 3 is also used to measure the decay D* 2 ^ —*■ 

The detector resolution in the reconstruction of this decay is (3.1 ± 0.1)MeV/c^. The 
decay is also allowed but suppressed due to the smaller phase space. 

The measured mass-difference distribution for K"*" combinations with a scaled 
energy greater than 0.25 is shown in Figs. ^-f. Signals for the Djj ^ (Fig. ^) 












Decay channel 

c-quark-enriched sample 

b-quark-enriched sample 

D+ ^ D*+K° 

21.4 ± 5.0stat± 0.2gyst 

16.1 ± 5.3., 

D+i ^ D*°K+ 

32.2 ± 9.3stat± l-Ssyst 

15.4 ± 7.8stat± l-Ssyst 

D: 2 + ^ D°K+ 

30.7±14.Wlil,„ 

14.6±14.2..„«|,^,, 


Table 5: Number of signal events in the c- and b-quark-enriched event sample. 


and D*^ —> decays (Fig. are visible, as expected, around 525 MeV/c^and 

710 MeV/c^, respectively. The shape of the combinatorial background is determined with 
the help of the wrong-sign charge combinations D° K“(Figs. |^d and f). The shape of the 
combinatorial background for the decay of the is parameterized with the function given 
in Section PTT| . The combinatorial background for the decay of the D*^ is described by a 
linear function of AM 3 . The numbers of signal events are determined with a simultaneous 
unbinned log-likelihood £t to the signal and background distributions. The signal is 
fitted with a Breit-Wigner function convoluted with the Gaussian detector resolution. 
In addition, reflections from the decays and D 2 ''' ^ are taken into 

account in the £t to the and distributions. This contribution 

is small due to the broad distribution of the reflection events and is determined from 
simulated events to be 4% ( 6 %) of the ^ D*°K+ (D*^ —> D°K+) signal events. 
The results of the £t to the AM 3 distributions are listed in Table The systematic 
uncertainties on the measured mass differences are obtained by varying, in the £t, the 


assumed decay width and the detector resolution, as described in Section |5T . 

From the measured mass differences, the mass of the and can be calculated 
by adding the appropriate nominal mass (D*’*', or D*’) and calculating the weighted 
average in the case of the D^. The resulting masses are m(D^)=2535.3 ± 0.7 MeV/c^ 
and m(D*^)=2568.6 ± 3.2 MeV/c^. The measured masses are in good agreement with 
the world averages given in Table 0. 


3.3 Measurement of D** decays in the c- and b-quark-enriched 
sample 

The measured mass-difference distributions for the c- and b-quark-enriched samples are 
shown in Figs. ^ and |^. The Dgi and D *2 mesons are reconstructed as described in the 
previous sections. The numbers of signal events are determined with an unbinned log- 
likelihood fit to the data as described in the previous section, with the exception that the 
mass differences are hxed using the Dsi and D *2 masses given in Table ^ The results are 
given in Table |^. 


4 Reconstruction efficiencies 

To determine the production rates in hadronic events, the reconstruction efficiencies 
are taken from the simulation. The efficiencies are calculated separately for D** mesons 
produced in charm fragmentation, b-hadron decays and from gluon splitting. 

In the flavour-independent selection the reconstruction efficiencies for D** mesons 
produced in charm fragmentation and in b-hadron decays are approximately the same. 
The efficiency for D** mesons from gluon splitting is lower due to their softer energy 
spectrum. The mean efficiencies are calculated from the ffavour-specihc efficiencies under 
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Figure 4: Reconstructed mass-difference distribution for the c-quark-enriched event 
sample. The full curve is the result of an unbinned log-likelihood £t and the dashed 
curve represents the combinatorial background. 


The 


the assumption that D** mesons are produced in equal amounts in charm fragmentation 
and in b-hadron decays. The fraction of D** mesons produced in gluon splitting is taken 
from the simulation, with the rate of gluons splitting into cc pairs reweighted to the 
measured value of /(g —>• cc )=0.0296 ± 0.0038 ||I5| 
average efficiency is calculated as = 0.47 ■ ec -|- 0.47 
efficiencies are given in Table p. 

The reconstruction efficiencies in the c- and the b-quark-enriched event samples for 
the different D** decays are listed in Table |^. 


and found to be about 6 
eb -l- 0.06 


eg. The calculated mean 


5 Production rates 

The production rates are calculated from the observed number of events, the total number 
of hadronic Z decays (A'"evt=4 151 890), the known branching ratios and the reconstruction 
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Figure 5: Reconstructed mass-difference distribution for the b-quark-enriched event 
sample. The full curve is the result of an unbinned log-likelihood £t to the data and 
the dashed curve represents the combinatorial background. 


efficiencies 

/(Z ^ D±) ■ Br(D+ ^ D*+K0) 
/(Z D±) ■ Br(D+ ^ D*0K+) 

/(Z ^ ) ■ Br(D:+ ^ D0K+) 


2 ■ Til 


iVevt • ei ■ Br(D*+ ^ D°7r+) ■ Br(K0 ^ tt+tt-) 
_ ^2 _ 

iVevt • €2 ■ Br(D0 ^ K-7r+) ’ 

_ ^ _ 

iVevt • 63 • Br(D0 ^ K-7r+) ’ 


( 1 ) 

( 2 ) 

(3) 


where ui, U 2 and denote the number of reconstructed events in the three considered 
D** channels and ei, 62 and 63 the corresponding branching-ratio-weighted sums of the 
reconstruction efficiencies. The branching ratios of the D*"*", Kg, D°, and those of all 
subsequent decays, are taken from [ 0 . 

In the case of the c- and b-quark-enriched samples, the rate is calculated per event 
hemisphere and not per Z-decay. In this case the number of hadronic events are multiplied 
by the partial width Rc and Rb, for which the Standard Model values Rc = 17.2% and 
















































Decay channel 

Mean efficiency 

Dg+ ^ D*+K0 (D° ^ K-7r+) 

(16.46 ±0.41)% 

Dg+i ^ D*+K0 (D° ^ kV+tt-) 

(11.70 ±0.51)% 

Dg+i ^ D*+K0 (D° ^ K-7r+7r-7r+) 

(9.17 ±0.22)% 

Dg+i ^ D*+K0 (D° ^ K-tt+ttO) 

(2.99 ±0.11)% 

Dg+i ^ D*°K+ 

(11.00 ±0.14)% 

D:2+ ^ D°K+ 

(10.73 ±0.18)% 


Table 6 : Mean reconstruction efficiencies of the decays ^ and 

0^2+ ^ D°K+. 


Decay channel 

c-quark-enriched 

b-quark-enriched 

c ^ Dg** 

b ^ Dg** 

c ^ Dg** 

b ^ Dg** 

Dg+i ^ D*+K° 

(5.87 ±0.30)% 

(0.88 ±0.14)% 

(0.34 ±0.04)% 

(4.08 ± 0.20)% 

Dg+i ^ D*°K+ 

(9.14 ±0.22)% 

(1.40 ±0.06)% 

(0.70 ±0.06)% 

(7.81 ±0.15)% 

Dg* 2 + ^ D°K+ 

(9.99 ±0.18)% 

(1.36 ±0.08)% 

(0.59 ±0.05)% 

(7.34 ±0.17)% 


Table 7: Reconstruction efficiencies for the decays and 

Dg^ —> in the c- and b-quark-enriched event samples. For the decay 

the reconstruction efficiencies are branching-ratio-weighted sums of the four different 
D° decay modes. The low efficiency for b —^ D** (c —>■ D**) in the c-quark (b-quark) 
enriched sample reflects the contamination of approximately 10 — 15% of the other flavour. 


Rb = 21.6% are used. The average number of D** expected from gluon splitting (1 — 2%) 
is subtracted from the number of events observed. Since the enriched samples are not 
pure, but have a contamination of roughly 10 % from other quark flavours, the production 
rates are determined simultaneously for the c- and b-quark-enriched samples. 

5.1 Systematic errors 

The main contribution to the systematic error originates from the uncertainty on the 
parameters held hxed in the £t to the mass distributions of the and D*^ candidates. 
In all hts the natural decay widths (r(D^) = 1.5 MeV and r(D*^) = 15 MeV) and the 
detector resolution are hxed. In order to determine the systematic error, the decay width 
of the (Dg^) is varied between 1 and 2 MeV(10 and 20 MeV). An uncertainty of ±10% 
is assumed for the detector resolution determined from the simulation. In the case of the 
c- and b-quark-enriched samples, the central values of the and masses are hxed to 
the world averages 1^ and are varied by ±0.6 MeV/c^for the and by ±2 MeV/c^for 
the Dg^. Each parameter is varied independently and the change in the number of signal 
events is taken as the systematic error from that particular source. The total systematic 
error listed in Tables ^ and is given by the quadratic sum of all the above contributions. 

The rate of gluons splitting to cc pairs is varied within the experimental uncertainties. 
In addition, the selection efficiency for D** mesons originating from this source is varied 
by 50% to account for the modeling of the g —> cc process. An error of 20% is assumed 






















Systematic errors in the decay ] 

3*+K0 


Error source 

flavour independent 

c-enriched 

b-enriched 

Mass, decay width and resolution 

+8-9 

+ 1 

+5 -11 

xe{s), XEih) variation 

± 1 

± 1 

± 1 

Kaon dE/dx selection 

± 3 

± 3 

± 3 

D** source (gluon, c and b) 

± 3 

± 1 

± 1 

Monte Carlo statistics 

± 3 

± 3 

± 4 

branching ratios 

± 5 

± 5 

± 5 

Total 

+11 - 12 

± 7 

+9 -13 

Systematic errors in the decay ] 

3*ok+ 


Error source 

flavour independent 

c-enriched 

b-enriched 

Mass, decay width and resolution 

± 8 

± 4 

± 12 

xe{c), a;£;(b) variation 

± 1 

± 1 

± 2 

Kaon dE/dx selection 

it 3 

± 3 

± 3 

D** reflections 

± 2 

± 4 

+3-2 

D** source (gluon, c and b) 

+ 3 

± 1 

± 1 

Monte Calo statistics 

± 1 

± 1 

± 1 

branching ratio 

± 2 

± 2 

± 2 

Total 

± 10 

± 7 

± 13 

Systematic errors in the decay > 

D 0 K+ 


Error source 

flavour independent 

c-enriched 

b-enriched 

Mass, decay width and resolution 

+7 -15 

+9 -17 

+18 -34 

xe{c), xe{/>) variation 

± 1 

± 1 

± 2 

Kaon dE/dx selection 

± 3 

± 3 

± 3 

D** reflections 

± 3 

± 3 

± 14 

D** source (gluon, c and b) 

± 3 

± 1 

± 1 

Monte Carlo statistics 

± 1 

± 1 

± 1 

branching ratio 

+ 2 

± 2 

± 2 

Total 

+9 -16 

+10 -18 

+23 -37 


Table 8 : Systematic error contributions in %. 























































































































on the equality of the D** production in c and b events, which is used in the calculation 
of the average efficiency (see Section The resulting error on the production rate is 1%. 

The fraction of reflection events from D** decays is taken from the simulation. The 
obtained ratio between reflection and signal events is varied by 50% and the corresponding 
change in the number of signal events taken as the systematic error. 

In the simulation the parameters of the Peterson fragmentation function are 
reweighted to reproduce the measured values of the mean scaled energy xg of c 
and b hadrons: xe{c) = 0.484 ± 0.008 and X£;(b) = 0.702 ± 0.008 0) respectively. The 
reweighted samples are used to determine the reconstruction efficiencies and the errors on 
the fragmentation parameters are propagated as well. 

The dE/dx selection efficiency determined from the simulation is tested by 
reconstructing D*"*" —> (K“ 7 r’'') 7 r’'' decays without using the dE/dx information of the 

K“. By applying subsequently the dE/dx cut used in the analysis, the efficiency of this 
cut can be determined in data and the simulation. The difference between data and 
the simulation is (—3.8 ± 2.9)%. The reconstruction efficiencies of the appropriate decay 
channels are corrected for the observed offset and the uncertainty is taken into account 
as a systematic error. 

Finally the errors of the reconstruction efficiency due to the hnite Monte Carlo 
statistics and the uncertainty on the used branching ratios are taken into account. 

The total systematic error is calculated by adding the different contributions in 
quadrature. The list of the systematic errors is given in Table p. The systematic errors of 
the two Dsi-modes have been assumed to be 100% correlated. 


5.2 Results 

From the measured numbers of produced and D*^ mesons and the reconstruction 
efficiencies, the production rates in the decay channels ^ and 

D*^ —> in hadronic events are 

/(Z ^ D±) ■ Br(D+ ^ D*+K0) = (0.22 ± 0.05,tat± 0.03,y,0% 

/(Z ^ D±) ■ Br(D+ ^ D*0K+) = (0.29 ± 0.08,tat± 0.03,y,0% 

/(Z ^ D:±) • Br(D:2+ ^ D0K+) = (0.37 ± 0.13,tatlli:g^,,,J% • 

The production rates in the c-quark-enriched sample are 

/(c ^ D±) ■ Br(D+ ^ D*+K0) = (0.35 ± 0.10stat± 0.03syst)% 

/(c ^ D±) ■ Br(D+ ^ D*0K+) = (0.59 ± 0.19stat± 0.04y,t)% 

/(c ^ D:±) ■ Br(D: 2 + ^ D0K+) = (0.51 ± 0.27,tatlg:gEsyst)% • 

The production rates in the b-quark-enriched sample are 
/(b ^ D±) . Br(D+ ^ D*+K») = (0.30 ± 0. 

12stat-0.04syst)^ 

/(b ^ D±) ■ Br(D+ ^ D*0K+) = (0.24 ± 0.15stat± 0.03,yst)% 

/(b D: 2 ±) ■ Br(D: 2 + ^ D0K+) = (0.25 ± 0. 

29stat-0.09syst)^ • 

The production rate of mesons in b-decays has a signihcance of less than one sigma, 
therefore a 95% conhdence limit is calculated by renormalizing the probability for the 
allowed region of the production rate and taking into account the measured central value: 

/(b ^ ■ Br(D: 2 + ^ D°K+) < 0.77% at 95% CL . 




Assuming that the decay width of the is saturated by the two measured decays and 
that the branching ratio Br(D*^ —> is 45%, which follows from the assumption 

that Br(D*^ ^ D K)=90% with 10% branching ratio into D* K [||, the total production 
rates in hadronic events are found to be 

/(Z ^ D±) = (0.52 ± 0.09stat± 0.06sy«t)% 

/(Z ^ D:2±) = (0.83±0.29stattHLyst)% ■ 

The production rates from c-quarks (c —> D**) are 

/(c ^ D±) = (0.94 ± 0.22stat± 0.07,y,t)% 

/(c ^ D: 2 ±) = (1.14±0.59,tattliiJ,y,J% . 

The production rates from b-quarks (b —> D**) are 

/(b ^ D±) = (0.55±0.19.„^g:g?,J,„)% 

/(b ^ d:±) = (0.57 ± 0.63.,..ig:gg,j„)% . 

The production rate of D*^ mesons in b-decays has a signihcance of less than one sigma 
and the corresponding upper limit is 

/(b ^ D: 2 ^) < 1.70% at 95% CL . 


6 Conclusions 


The production rates of the D** mesons and have been measured by reconstructing 

the decays D+i ^ D*+K0 , D+i ^ D*^K+ and 0 ^ 2 + ^ D0K+. 

In the two measured decay modes an enhancement of the hnal state is 

expected due to the higher Q-value of the decay, since isospin invariance requires the 
matrix elements of the two decays to be the same. The expected ratio of the two branching 
ratios is given by P] 


Br(D.+i ^ D*°K+) _ /q(KU\“+' 

Br(D+ ^ D*+K0) q(K0 )) ~ ’ 

where q(K+) (q(K°)) is the momentum of the K’''(K°) in the rest frame of the and 
L=2 for a pure D-wave decay. 

The measured ratio of the two production rates is 


Br(D+i ^ D*°K+) 
Br(D+i ^ D*+K°) 


1.32 ±0.47 


stat 


± 0.23syst. 


The smaller value compared with the expectation could be due to an admixture of a 
S-wave decay from one of the broad D** states (see Table |I|). However, no signihcant 
discrepancy with the expected production rate ratio for a pure D-wave decay is observed. 
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